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We study the effect of quenched spatial disorder on the steady states of driven systems of in- 
teracting particles. Two sorts of models are studied: disordered drop-push processes and their 
generalizations, and the disordered asymmetric simple exclusion process. We write down the exact 
steady-state measure, and consequently a number of physical quantities explicitly, for the drop-push 
dynamics in any dimensions for arbitrary disorder. We find that three qualitatively different regimes 
of behaviour are possible in 1-d disordered driven systems. In the Vanishing-Current regime, the 
steady-state current approaches zero in the thermodynamic limit. A system with a non-zero cur- 
rent can either be in the Homogeneous regime, chracterized by a single macroscopic density, or the 
Segregated-Density regime, with macroscopic regions of different densities. We comment on certain 
important constraints to be taken care of in any field theory of disordered systems. 
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I. INTRODUCTION 

It is known that quenched disorder can strongly af- 
fect the large-scale, long-time behaviour of nonequilib- 
rium driven systems with interacting constituents. The 
interplay of disorder, interactions and drive opens up the 
possibility of new regimes of complex and interesting be- 
haviour arising in these systems [Q. In the theoretical 
effort to delineate and explore regimes of new behaviour, 
an important role is played by simple models which cap- 
ture some features of more complex physical systems. In 
this paper, we study disordered driven diffusive systems 
by analysing stochastically evolving lattice gas models, 
with quenched disordered hopping rates ||| . 

Driven diffusive systems in the absence of disorder 
have been studied extensively and are reviewed in ||. 
Also, systems with disorder and drive but no interac- 
tions between particles are well studied and understood 
||. But there have been only sporadic studies of dis- 
ordered driven diffusive systems of interacting particles. 
It has been argued that strong enough random site di- 
lution can substantially affect the transport properties 
of particles with hard-core interactions, and can make 
the system respond nonmonotonically to the driving field 
HP)- O n the other hand, a low concentration of blocked 
sites was found numerically not to affect the critical be- 
haviour of a driven lattice gas with additional attractive 
inter-particle interactions M. Finally, a driven lattice 
gas with a quenched noise distribution was studied us- 
ing field-theoretic techniques in ||, but the connection 
of this study with particle-conserving disordered lattice 
gas models is not clear. 

In this paper, we study disordered lattice gas models 
with a view towards identifying different sorts of generic 
behaviour that can arise on large scales as a consequence 
of disorder. The only interaction included is the hard- 
core constraint which limits the allowed occupancy of 



each site. Our results pertain mostly, but not exclusively, 
to one dimension. In the remainder of this Introduction, 
we discuss the different types of behaviour displayed by 
the lattice gas models under study. 

We find three distinct regimes in disordered driven dif- 
fusive systems in one dimension: 

In the Homogeneous regime, the state of the system is 
characterized by a single density and a nonzero current. 
Quenched disorder induces variations of the density on 
the microscopic scale, of the order of a few lattice spac- 
ings. However, the system has a macroscopically homo- 
geneous density. In the thermodynamic limit, the current 
approaches a finite value. 

In the Segregated-Density regime, the state of the sys- 
tem is characterized by two distinct values of density, and 
a nonzero current. Besides microscopic-scale variations of 
the density, there are macroscopic regions with differing 
high and low densities. The state is thus characterized 
by phase separation of the density, and a spatially con- 
stant time-averaged current which remains finite in the 
thermodynamic limit. 

In the Vanishing- Current regime, the state of the sys- 
tem is characterized by two distinct values of the den- 
sity, and an essentially zero current. The hallmark of 
this regime is that the current decreases as the system 
size increases, and vanishes in the thermodynamic limit. 
This is a consequence of rare but rate- limiting backbends, 
or stretches of bonds which disfavour the forward flow of 
current. The density is inhomogeneous on a macroscopic 
scale. 

The density profiles in typical states in each of the 
three regimes are depicted in Figure f, while Figure 2 
shows the variation of the current with system size in 
the three cases. 
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FIG. 1. Representative steady-state density profiles for the 
(a) Homogeneous (b) Segregated-Density and (c) Vanish- 
ing-Current regimes in the Disordered Asymmetric Exclusion 
Process (DASEP). 
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FIG. 2. Variation of the steady-state current with the sys- 
tem size for the three DASEP regimes of Fig. 1: (a) Homo- 
geneous (circles), (b) Segregated-Density (triangles) and (c) 
Vanishing-Current (squares). In (a) and (b), the current ap- 
proaches a finite value in the thermodynamic limit whereas 
in (c) the current vanishes as a power of the system size. The 
dashed line corresponds to J = 0.125 which is the limiting 
value of the current in the regime (b) for the chosen values of 
the parameters. 

Examples of these behaviours are discussed in this 
paper for two types of lattice-gas models, namely the 
disordered drop-push process (DDPP) and the disor- 
dered asymmetric simple exclusion process (DASEP). 
The models are defined in detail in Sections || and III 
respectively, but for the purpose of discussion here, it suf- 
fices to note that the models are similar in that there is a 
maximum occupancy of each site in both, and are differ- 
ent in the dynamical moves — attempted nearest neigh- 
bour jumps in the DASEP, and slightly longer-ranged 
leapfrogging moves in the DDPP. 



The absence of detailed balance, together with the 
breaking of translational invariance, in disordered off- 
equilibrium systems makes the characterization of even 
the stationary state difficult in general. It is shown that 
the steady state of the disordered drop-push process can 
be found explicitly - the first such instance we are aware 
of, in a system with disorder, interactions and drive ||. 
This determination - which is based on the condition of 
pairwise balance [^0) - shows that a product measure 
form is valid in all dimensions. The form reflects the 
microscopic inhomogeneities coming from the underlying 
disorder, and results in a macroscopically homogeneous 
state. 

For the disordered asymmetric exclusion process, the 
steady state measure is not analytically characterizablc, 
and we study the problem within a site-wise inhomo- 
geneous mean-field theory and by numerical simulation. 
The result depends crucially on whether or not the sys- 
tem has backbends, which are stretches of the lattice 
where the local bias is against the particle flow. In the 
no-backbend case, when the average particle density is 
sufficiently away from 1/2, the spatial profile of the den- 
sity has microscopic shocks, but is uniform on macro- 
scopic scales (Fig. |I|a). However, in a finite region around 
half-filling, disorder induces phase separation into macro- 
scopic regions of high and low density (Fig. |l]b) . We give 
approximate arguments to understand the origin and na- 
ture of this phase separation, and to obtain the form of 
the phase diagram in the current-density plane. This sort 
of behaviour has also been seen earlier in a model with 
a single weak bond |ll| . We argue that disorder- induced 
phase separation is a generic feature of systems in which 
the current J versus density p shows a maximum at some 
intermediate density, in the absence of disorder. 

In the version of the DASEP in which the easy di- 
rection of hopping is itself a quenched random variable, 
the model represents a system of hard-core particles in a 
random potential with an overall downward tilt, but with 
backbends of arbitrary length. Long backbends severely 
limit the maximum current that can flow through the sys- 
tem, and in fact the current decreases to zero as the sys- 
tem size increases (Fig. ^]); the system is in the vanishing- 
current regime. 

Although our emphasis in this paper is on the analysis 
of lattice models, we comment briefly on certain con- 
straints that are important in a continuum description. 
Such a description is expected to be valid for the large- 
scale, long-time behaviour, and is based on stochastic dif- 
ferential equations involving appropriate coarse-grained 
variables. It is argued that quenched randomness is man- 
ifest in random multiplicative coefficients in a gradient 
expansion. Conservation of particle number - which im- 
plies spatial constancy of the current in the steady state 
- imposes strong constraints on these terms. 

In one dimension, using a well known mapping [ jl2| , 
the particle models are equivalent to stochastic growth 
models of a 1-d interface moving in a 2-d medium. The 
interface moves with a speed proportional to the cur- 
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rent in the particle model. The disordered jump rates 
now become local growth rates which are disordered in a 
columnar fashion for the moving interface [|l3| . The three 
principal regimes of behaviour discussed above for the 
particle models translate into distinct regimes for inter- 
face motion, namely (i) a moving interface with normal 
roughness, (ii) a moving interface with large segments 
with different mean slopes, and (iii) an interface with 
different-slope segments, which is stationary in the ther- 
modynamic limit. 

The paper is organized as follows. In Section || wc 
define and discuss the steady state properties of the dis- 
ordered drop-push process in arbitrary dimensionality. 
The disordered asymmetric exclusion process with only 
forward-easy-direction of hopping, but quenched random 
rates, is discussed in Section III; the case in which there 



are some backward-easy-direction bonds is discussed in 
Section [V. In Section [V] we discuss the constraints on 
a continuum description, while Section VI discusses the 



implications of our results for models interfa ce g rowth 
in the presence of columnar disorder. Section VII is the 
conclusion. 



II. DISORDERED DROP-PUSH PROCESS : 
DDPP 

The drop-push process was initially introduced in 
P,^0| as a model of activated flow involving transport 
through a series of traps of equal depths. The dynam- 
ics consists of activated hops together with a cascade of 
overflows following each move. The disordered version 
of the model may be considered as a discrete model of 
activated fluid flow down an inclined rugged slope with 
lakes of varying depths; see Figs. |L|| This is similar 
to above-threshold behaviour of the model considered in 
||l4|| . In this section we show that the steady state and 
current can be found exactly in all dimensions for the 
DDPP and its generalizations. 



FIG. 4. A Disordered Drop-Push Process (DDPP) config- 
uration and move in d = 1. 



A. The model 

The model in d-dimensions is defined on a hypercubic 
lattice with periodic boundary conditions along all the d 
axes (with unit vectors {e v \v = 1 • • • d}). At each site r is 
a well which can hold at most l r particles (Figs. with 
Z r 's chosen independently from some probability distribu- 
tion P(l). The configuration C of the system is specified 
by specifying the set of occupation numbers {n r } with 
(0 < n r < l r , Vr). Further, with each site r is assigned a 
set {e(n r \l r ); n r = 1, • • • , l r } of positive random numbers 
chosen from some given distribution Jl5[ ] . The dynam- 
ics is stochastic. In a time interval dt, with a proba- 
bility p± v e(n r \l r )dt, the topmost particle in the well r 
hops out, and drops into well r ± e„, i.e. into the adja- 
cent well in the ±^th direction. Here {p± v ] v = 1, .., d} 
are a set of site-independent positive numbers satisfying 

J2t=i(P» + V-v) = 1- Now, if well r ± e„ is already full, 
then the particle gets pushed further preserving the di- 
rection of the initial jump to the next site and so on. 
The cascade of transfers terminates once a partially full 
well is encountered. Note that here the set of jump-rates 
{e(n r |Z T .)} are site-dependent as well as functions of the 
occupation numbers. These rates, together with the well- 
depths {l r }, constitute the quenched random variables in 
the model. The set of probabilities {p±v} determines the 
direction of the global bias E_^_J2t=i(.Pv ~ P-v)&v and, 
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as will be shown in Section [II B[ also the direction and 
magnitude of the steady state current in the model. How- 
ever, they do not enter the expression for the normalized 
invariant measure. 




FIG. 3. Schematic diagram of water flowing down a rugged 
hill-side. Water from a lake higher up cascades downhill, un- 
der the action of gravity, until it finds a partially filled lake. 
The unequal capacities of the lakes are the quenched variables 
in the system. 
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FIG. 5. The DDPP model in d = 2. The model can be 
generalized to d > 2 (see text). The rates e(n r \l r ) depend on 
the well depth l r as well as the occupation number n r . 
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Though all the results we will discuss holds for any ar- 
bitrary choice of the e's, in a physical system they should 
be determined from the details of the trapping mecha- 
nisms etc., e.g. they may be taken to be of the Kramers 
form e(n r \l r ) oc exp\—g(l r — n r )] for situations where the 
jumps are activated |lq|. 



B. Invariant measure 

The time evolution of the probability V(C) for the sys- 
tem to be in configuration C is given by the master equa- 
tion iIItIi 



j f v{c) = w{ £" -» cyp{c") - w ( c -* c ')^(c). 

(1) 



Here the Ws are the transition matrix elements iden- 
tified with the rates e's defined in the model, e.g. if 
the transition C — > C involves moving the topmost of 
the n r particles at r to r' along the vth direction, then 
W(C — > C) = p v e r (n r \l r ). The steady-state or the in- 
variant measure of the dynamics is the set of time in- 
dependent weights {(J-(C)} satisfying ([!]) above. Hence 
the problem of finding the invariant measure reduces to 
that of finding a set of positive weights {/x(C)} such that 
the total incoming flux into any configuration C (the first 
sum in (0)) equals the total flux out of C (the second 
sum in (□})■ The uniqueness of the invariant measure is 
ensured by the connectedness property of the VF-matrix, 
i.e. every configuration can be reached from any other 
by a sequence of transitions jlTj . 

We claim that the (unnormalized) measure of configu- 
ration C({n r }) in the steady state has the product form 

/i(C) = J[u r (n r ). (2) 

r 

Here u r are the single-site weights defined as 



u r (n r ) 



(1) • • -T r (n r ) 



if n r = 

if < n r < L 



(3) 



where T r (n r ) — eo/e(n r \l r ) with eo being a microscopic 
rate. 

To show that (|2|) is indeed the invariant measure for 
the DDPP we show that it is possible to associate con- 
figuration C" in one-to-one correspondence with every C 
obtained from C by an elementary transition such that 



W(C -> C')n(C) = W{C" -> C)^(C"). 



(4) 



The above is the condition of pairwise balance fl(|| which 
ensures that the terms in the two sums on the right hand 
side of (|l]) cancel in pairs. Pairwise balance has been 
used earlier to find steady states of translationally in- 
variant systems |lO|,[L8|]. We now see that it can be used 



effectively to deduce the steady state of a disordered sys- 
tem as well. 

Suppose the transition C — * C involves hopping a par- 



ticle at site r to a site r' 



Ar'e„ with all wells in 



between along the vth axis full (Fig. |6|). Also suppose the 
well r" = r — Ar"e„ is not full but all wells between r" 
and r are full. The configuration C" is constructed such 
that it is identical to C at all sites except at the sites r" 
and r, at which n"„ = n r n + 1, n'J. — n r — 1. With C" 
thus defined, the pairwise balance condition (||) 



t(n r \l r ) Y[ u r (n r ) — e(n"„ \l r ») u r {n'^ 



(5) 



reduces to 



e(n r \l r )u r "(n r ")u r (n r ) — e(n"„\l r ")u r "(n"„)u r (n"), (6) 

since C and C" differ only at the sites r" and r. As can be 
explicitly checked, condition (||) is satisfied in view of the 
form of the weights u r (n r ). The prescription above en- 
sures one-to-one correspondence between configurations 
C" and C — > C transitions. 



C" 



~1M' 



C 



FIG. 6. Construction of configurations satisfying the pair- 
wise balance condition in d — 1. 

In the limit of large L we use the grand canonical de- 
scription and write the normalized probability of config- 
uration C in the steady state as 



V(C) 



//(C) 



UrZr 



where the site generating functions are given by 



n,.=0 



u r (n r )z nr 



(7) 



(8) 



and N p is the number of particles in the configuration. 
Here z is the fugacity and will be shown to be directly 
related to the steady-state current. The mean particle 
density at site r is given as (n r ) = zdlnZ r /dz Thus 
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the fugacity z is related to the mean particle density p of 
the system through 



z sr-^ dlnZ r 



dz 



(9) 



The steady state is characterized by a spatially uniform 
z but inhomogeneous site densities. 



C. Steady-state current 

For simplicity of presentation, we first derive a closed 
form expression for the steady-state current Jo for the 
1-d fully asymmetric model. Then we use these results 
to write down the current for the general ci-dimensional 
case. 

(i) d — 1, with forward jumps: In this case, for nota- 
tional convenience, we replace the lattice site index r by 
a single integer index i. Further, we allow jumps only in 
one direction (y = 1 : pi = p,p-i — 0). For an infinite 
system we write the current Jj^+i in the bond (i, i + 1) 
as 

h 

Ji,i+i = E E P e ( n i\h) p A n j) II Pk ^ 

j<i rij=l j<k<i 



where 



2J pe(ni\li)Pi(rii) 



Pi{rii) = Ui(ni)z ni /Zi 



(10) 



(11) 



is the probability that site i has occupancy rij. Note that 
the product form of the measure (^) has been exploited 
to write the above expression with decoupled joint prob- 
abilities. 

It is cumbersome, and not very instructive, to perform 
the sum in ( [To| ) to obtain the current in a closed form. 
Instead we use the spatial constancy of the current in the 
steady state to arrive at the result. We note that in the 
above expression the first term represents contribution of 
jumps originating from site j to the left of site i with all 
sites in between full. The second sum represents jumps 
from the site i itself contributing to J^j+i. The first sum 
may be rewritten as 

+ Y,n7-i=l e ( n i-l\h-l)Pi-l(rii-l) ■ 

The quantity within square brackets [•] is immediately 
recognised to be Jj-i,, by comparing it with (|io|). Physi- 
cally, this means that all jumps contributing to Ji-i,,:, for 
which site i is completely full, also contribute to J^j+i. 
Hence we have the recursion 



Ji 



Pi(k)Ji-i,i + V pe(ni\li)Pi(ni) (12) 



relating Jj.i+i to Jj-i^. Now, since {Pi(rii)} are the 
steady-state site probabilities, and since in the steady 
state all the bond currents must be equal (i.e. Jij+i = 
Ji-i,i = = Jo), from (|l2| ) we obtain 



■h 



T,n, = l e ( n i\ l i) P i( n i) 
l-Pi{k) 



pe z 



(13) 



where (g) and (^l|) have been used in the second step 
above. 

Note that the steady-state current does not depend 
upon the detailed spatial arrangement of wells. It is only 
a function, through the fugacity z, of the density and the 
total number of wells of different types in a particular 
realization of disorder. 

(ii) d = 1, with jumps in both directions: We can 
write the current J^i+i in bond (i, i + 1) as the difference 
between the current J[ i+ i due to the rightward jumps 
and the current due to the leftward jumps. This 
can be done since in each cascade the direction of the 
initial jump is preserved. Now, we use the result for the 
fully asymmetric case above to each of these currents 
separately to obtain J^j+i = Jo = {pi — P-i)€qz 

(iii) d > 1: To generalize the above results to d > 1 
we note that for DDPP in d > 1, the invariant measure 
(Q) is the same if we single out a particular direction, say 
v , and allow jumps only along that direction. Together 
with the direction preservation of individual jumps, this 
allows us to write the expression for the current in any 
dimension: 



Jo 



P-v)e v 



e z = e zE 



(14) 



where E = ^2 u —i(p v — P-v)e v is the external drive. As 
in the d = 1 case, the magnitude and the direction of the 
steady-state current does not depend upon the detailed 
arrangements of the wells. 



D. Static two-point correlation functions 

Because of the product form of the measure, the con- 
nected part of the equal-time density-density correlation 
function 



G r (Ar) = (n r n r+ Ar) - (n r ){n r +Ar) 



(15) 



vanishes identically for Ar ^ 0. Consequently, the fluc- 
tuation of the number of particles in r consecutive sites 
along a straight line can be computed exactly: 



i+r 
j=i+l 



<»i» 



i-\-r 



- e (4) **>■ 

3 =i+l 



(16) 
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The second step follows from the product form of the 
measure. 

For d = 1, a standard mapping discussed in Section VI 
introduces height variables defined by 



hi =^T f 2({nj) -rij). 

j<i 



(17) 



Evidently, r?(r) is the equal-time height-height corre- 
lation ((hi+ r — hi) 2 ). Averaging over the disorder dis- 
tribution gives r 2 (r) <~ r implying that the 'roughness' 
exponent a (defined by T 2 (r) ~ r 2a ) is 1/2. 



E. Two-rate DDPP model: Explicit results 



Let us consider a drop-push model where the maximum 
occupancy of each site is restricted to one, i.e. l r = l,Vr, 
but the hopping rates e r (l) are disordered, and chosen 
independently from the binary distribution 

Prob(x = e a ) = l-f, Prob{x = e b ) = /. (18) 

This model has the essential ingredients of disorder 
present in the original DDPP, yet it is simple enough that 
explicit, closed form relations between the mean density 
p and fugacity z, and hence the steady-state current Jo, 
can be written down. 

Let us denote by Z a and the site generating func- 
tions for the a and the b sites respectively. Using 
and (H), these are given by Z a = 1 + eoz/e a and Zt, = 
1 + eqz/ei,. Now, since the fractions of a and b sites are 
1 — / and / respectively, (0) reduces to 



Generalized Disordered Drop-push Process: 
GDDP 



We may consider a generalized version of the drop- 
push process in which, in addition to the particle moves, 
independent hole moves are also allowed. For simplic- 
ity we restrict ourselves to the generalized version of the 
single occupancy DDPP introduced above. This general- 
ized model may be regarded as the disordered lattice gas 
analo gue of the Toom interface dynamics in the low-noise 
limit 



see Section VI 



E Z 



Eb + e Z 



(19) 



The techniques developed for 
the DDPP may be used to obtain the exact steady-state 
measure and other quantities such as current and static 
correlations provided a certain condition [([22]) below] is 
met. 

The model in d dimensions is defined on a hypercubic 
lattice with periodic boundary conditions along all the d 
axes (with unit vectors \e v \v = 1 • • -d}), similar to the 
DDPP. Each site r of the lattice can hold either a parti- 
cle (n r — 1) or a hole (n r = 0). The configuration C of 
the system is specified by specifying the occupation num- 
ber of each well {n r } with (n r £ {0, l},Vr). Further, to 
each site r is assigned a pair of positive random numbers 
(a r ,0 r ) chosen from some distribution. The dynamics 
is stochastic and is very similar to that for the DDPP 
dynamics: in a time interval dt, a particle at site r is 
exchanged with the closest hole in the ±fth direction 
with a probability p± u a r dt (Fig. |?]). For identical parti- 
cles this move is equivalent to a cascade of particle moves 
terminating at the first vacant site as in the drop-push 
dynamics. Likewise, in interval dt, a hole at site r is ex- 
changed with the closest particle along the vth direction 
with probability q u j3 r dt. This can be looked upon of as a 
cascade of hole-moves analogous to the cascade of parti- 
cle moves. Here, as in the DDPP, the p v 's and g„'s are all 
non-negative and satisfy X^=i( a ^ + a -v) — 1; a = p,q- 
Further we chose a's and j3's such that 



This can be easily inverted to obtain z as a function of 
p, e.g. for / = 1/2 and E a = eq = Eb/q we obtain 



a r j3 r = K, 



(22) 



r x _ y/(l-g) 2 (l/2-p) 2 + g -(l + g )(l/2-p) 
(P) ~ 2(1 -p) [M) 

Since z(p) is known, the steady-state current is trivially 
obtained from (|l4|). 

Finally, the correlation function Tf(r) of (pi|), upon 
disorder averaging, may be written as 



( ^ ) [(1 -/)li)Z„+/h.Z,,] 



qz 



(I + 2) 2 (q + z) 2 
where z is given by (|20|) above. 



(21) 



o — o 




-o — o 



FIG. 7. Generalized Disordered Drop-push Process 
(GDDP) configuration and moves in d = 2. 
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where K is a constant independent of r. As we will 
see below, this particular choice of the jump-rates allows 
the exact determination of the invariant measure of the 
model. Physically this choice is quite reasonable since 
it implies that the sites which act as traps for particles 
(low a r ) are more transparent to holes (high (3 r ) and vice 
versa. As pointed out earlier the non-disordered version 
of this model, i.e. {a r = a, /3 r = /3, Vr} in l-d, is the lat- 
tice gas equivalent of the low-noise driven Toom interface 
dynamics 

The master equation (|l|) governing the time evolution 
of the system now includes terms corresponding to hole 
moves as well. Since each micro-step involves cither only 
particle moves or hole moves, we use the principle of pair- 
wise balance for the particle moves and hole moves sep- 
arately. We work in the grand canonical picture in the 
thermodynamic limit. 

If only particle moves were allowed the invariant mea- 
sure would be given by 



(C) = JJu r (n r ) 



f^pcle 



The single site weights u r are given by 



u r {n r ) 



1 

e /a r 



if n r = 
if rir = 1 



(23) 



(24) 



Introducing the fugacity z and site generating functions 
Z r = 1 + zu r (l) we can write the normalized single 
site probabilities as P r pde (0) = 1/Z r and PP cle {l) = 
zu r {\)/Z. r . 

Similarly, with only the hole moves, the invariant mea- 
sure has the product form 



Hhole(C) = Y\_ V r( n r), 



(25) 



where n r = 1 (0) refers to the presence (absence) of a 
hole. The single site weights v r are given by 



v r (n r ) = 



1 



eo//3r 



if n r 
if n r 



(26) 



Introducing the fugacity y for holes and site generat- 
ing functions Y r = 1 + yv r (l) we can write the nor- 
malized single-site probabilities as: P^ ole (0) = l/Y r and 
P r hole (l) = yv r (l)/Y r . 

Now, since each site is occupied either by a parti- 
cle or a hole, we must have PP cle (0) — P^ ole (l) and 
pP cle (l) = Pf° le (0). Using the detailed forms of PP cle, s 
and P^ ole 's, we arrive at the condition ( p2] ) with K = 
e$yz. If this condition is satisfied then the invariant mea- 
sure for GDDP is given by either ( p3f ) or ( |25| ) , since both 
are equivalent. 

In a similar manner as for the DDPP the current due 
to the particle moves and hole moves may be computed. 
The total particle current due to both types of moves, in 
d dimensions, is given by 



Jo = epz}^ -y- v )e v - e 2/^(g„ 



q-u)e v (27) 



As for the DDPP, static density-density correlations 
in the steady state vanish identically on account of the 
product form of the steady-state measure. In d = 1 the 
height-height correlation is given by ( |l6| ) . 



III. DISORDERED ASYMMETRIC SIMPLE 
EXCLUSION PROCESS: DASEP 

The asymmetric simple exclusion process (ASEP) is a 
prototype model for studying noncquilibrium phenomena 
in the context of lattice gases When discussing 

the effect of quenched disorder, it is important to dis- 
tinguish between cases in which (a) the easy direction 
of hopping in each bond is the same but hopping rates 
are random, and (b) the easy direction is itself a random 
variable. The latter case is studied in Section IV. In this 



section, we consider a 1-d system with disorder of type 
(a) and show that quenched disorder can induce macro- 
scopic phase separation. Using a variety of arguments we 
sketch the phase coexistence curve in the current {Jo) - 
mean density (p) plane. This agrees qualitatively with 
the results obtained from the Monte Carlo (MC) simula- 
tions. 



A. The model 

In one dimension, we define the disordered asymmetric 
simple exclusion process on a ring of L sites. Each site 
can hold either 1 or particle. To each bond + 1) 
of the lattice is assigned a quenched random rate a^j+i 
chosen independently from some probability distribution 
Prob(a). The dynamics is stochastic: in a time interval 
dt a particle at site i attempts to hop, with probability 
pai^+idt, to site i + The move is completed if and only 
if site i + 1 is unoccupied (see Fig. pf). 



a 



i,i+l 



- - o 



-o 



■o 



i i+1 

FIG. 8. The disordered fully asymmetric simple exclusion 
process (DASEP) and moves in d = 1. 

For the model defined above no analytically exact char- 
acterization of the steady-state measure could be ob- 
tained. A simpler model in which there is only one de- 
fect bond has been studied in detail by Janowsky and 
Lebowitz |Tl|] , but in this case too the exact steady-state 
measure is not known. We use Monte Carlo simulations 
and a mean-field approximation to demonstrate some 
striking effects of quenched disorder. 
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B. Current-density plot and density profile in steady 

state 

Figure ^| shows the steady-state current Jo vs mean 
density p plot, obtained from MC simulations, for a sys- 
tem of size L — 8000 and the rates a chosen from the 
binary distribution 



Prob(a = r) = f, Prob(a = !) = ! — /. 



(28) 
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B, which is a function of r and /, is obtained in subsec- 
tion III E below. 

We studied the steady state density profiles charac- 
terized by the set of site densities {pi = (rti)} in both 
regimes A and B, using MC simulations. We found that 
in regime A the system is homogeneous on a macroscopic 
scale, while in regime B it shows macroscopic density seg- 
regation. Figure [h] shows the steady state density pro- 
files for three representative mean densities — one from 
regime A and two from regime B. Evidently there is a 
large qualitative difference between the profiles in the 
two regimes. In regime A (Fig. fl0|a), there are density 
variations (shocks) only over microscopic scales; coarse- 
graining over a few lattice spacings leads to a spatially 
uniform density. In contrast to this, in the profiles corre- 
sponding to regime B (Fig. |lO|b), there are density inho- 
mogeneities over length scales comparable to the system 
size L, in addition to the shocks on a microscopic scale. 
This segregation into high and low density phases, with 
large shocks separating them, is reminiscent of phase sep- 
aration, and occurs over the full range of mean particle 
densities corresponding to regime B. A qualitatively sim- 
ilar phenomenon has been found in a system with one 
defect bond, studied in IllJl. 



FIG. 9. The current-density plot for the DASEP for a 
given realization of disorder for a system of size L = 8000. 
The hopping rates are chosen from the distribution ( p8[ ) with 
r = / = 1/2. The filled circles are the MC results and the 
dashed line is the mean field curve. The solid line represents 
the J — p curve for the Fully Segregated Model (FSM) for the 
same values of the parameters. 

Here / is the fraction of weak bonds, and r is a mea- 
sure of the strength of the weak bonds. We used the 
values r — 1/2, / =1/2 in our numerical work. For a 
specified mean density p, a random initial distribution of 
N p — pL particles on L sites is chosen and the system 
is allowed to settle into a steady state by evolving it for 
a sufficiently large number of MC steps. Then the cur- 
rent across each bond is obtained by counting the total 
number of jumps across that bond, over a large number 
of MC steps. An average of all the bond currents thus 
computed is taken as Jo, as currents across all bonds are 
equal in the steady state. Jo is a symmetric function of 
density around p = 1/2 as a result of a certain symmetry 
with respect to particle-hole interchange (see Appendix 
|a|). As may be expected, the current for the disordered 
system lies between the corresponding values of the two 
pure reference systems with r = 1 and r — 1/2 on all 
the bonds. The more striking qualitative difference be- 
tween the disordered and pure systems is the appearance 
of a plateau (Regime B in Fig. ^) for a range of densities 
— 1/2 1 < A. In this regime, the current is independent 
of the mean density and equals the maximum allowed 
current in the system. The approximate size A of regime 



0.5 

1 









FIG. 10. Density profiles for the DASEP for a system of 
size L = 8000 for a given realization of disorder at three fill- 
ings (a) p = 0.8, (b) p = 0.6, and (c) p = 0.5. In d,e and f 
are shown the blow-ups of the regions enclosed in the dashed 
boxes in a,b and c respectively. Circles are MC profiles and 
the continuous lines are mean-field results. 

The number of different large-scale regions of high and 
low density shows fluctuations from one realization of dis- 
order to another. As the size of the system is increased, 
we monitored the mean number of these regions, and 
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found that it is nearly constant, or perhaps increases very 
mildly - certainly, much less fast than linearly in the size 
of the system. This implies that the characteristic length 
scale of density segregation increases indefinitely in the 
thermodynamic limit. 



C. Mean- field approximation 

We now turn to a mean-field approximation, which as- 
sumes no correlations between site densities in the steady 
state. We will see that it captures most of the steady- 
state features found in the MC simulations above, includ- 
ing the steady-state density profile. 

The time-averaged steady-state current J%,%+\ in the 
bond is given by J M+ i = a M+ i(n; (l-n i+ i)). In 

view of the mean-field approximation (riirij) = (ni)(nj) 
this reduces to 



J, 



ati,i+i(rii) (1 - n i+ i) 



(29) 



To compute the steady-state current Jo as a function of 
the mean particle density p for a given realization of dis- 
order, we use the two iteration schemes described below, 
which yield equivalent results. 

(i) Constant current iteration scheme. For a given sys- 
tem of size L and for a fixed value Jo for the current, we 
iterate the following set of equations 



Pi+i = 1 - Jo/{oti,i+ipi) i=l,..,L 



(30) 



around the chain starting with, say, some value p\ (pe- 
riodic boundary conditions imply pi+L = Pi)- If Jo is 
less than a certain value J^axi which is the maximum 
current supported by the system within the mean-field 
approximation, then the iteration scheme converges, i.e. 
we get all the site densities in the physically acceptable 
range [0, 1]. The average of these site densities gives the 
mean density of the system corresponding to the station- 
ary current Jo. There are in general two values of the 
mean particle density corresponding to an allowed value 
of Jo and the iteration scheme converges to one or the 
other depending on the initial value of the density p\ . 

However, in this scheme we find that the number 
of iterations required for convergence increases without 
bound, as the trial value Jo gets closer to J^ax fr° m 
below. This divergence of the iteration scheme is pre- 
sumably due to the existence of the plateau in the J vs p 
curve, i.e. there exist many values of p for Jo very close 
to J^fx- Hence to obtain the density profile for p in 
the density segregated regime we resort to the constant 
density iteration scheme described below. 

(ii) Constant density iteration scheme. In this scheme 
we begin by assigning site densities {0 < Pi(0) < 1} to the 
lattice sites subjected to the constraint j^J2iPi(0) = P- 
The site densities are updated in parallel according to: 

Pi (t + 1) = Pi {t) + J, _!,*(*) - Ji,i+i(t); i = 1, ..,L 

(31) 



where J hi+ i(t) = a iii+1 pi(t)[l - p l+ i{t)] in view of 

We refer to this as the constant density scheme, since 
in each iteration the total density remains unchanged, i.e 
Si Pi(t + 1) = J2i Pi(t)- After a sufficient number of iter- 
ations, which depends upon the starting mean density p, 
the set of site densities converge to a set of numbers {pi} 
and the current on each bond converges to the steady 
state current Jo- 

The steady-state density profiles and the Jo vs p plot 
(0 < p < 1) for a given configuration of disorder obtained 
using these schemes is shown in Figs. ^ and ^| respec- 
tively. It is evident that the mean-field approximation 
( p9| ) reproduces quite well not only the J — p relation- 
ship, but also the density profile, including the locations 
of shocks, though not the shapes of individual shocks. 



D. Qualitative explanation of phase separation 



Although the mean-field approximation of the previous 
subsection successfully reproduces many features in the 
steady state, it does not yield a simple understanding of 
the phase separation (Fig. |l0|) or the plateau in the J-p 
curve (Fig. ^|) in terms of the macroscopic parameters of 
the model. We conjecture that underlying the behaviour 
of the DASEP in different regimes is a Maximum Current 
principle: For a given mean density, the system settles 
into a steady-state which maximizes the stationary cur- 
rent. Such a maximum current principle has been used 
to describe phase separation in the asymmetric simple 
exclusion process with open boundary conditions in [£3| . 

To use the maximum-current principle to have a qual- 
itative understanding of the phase separation in DASEP, 
let us assume that the density in each stretch of like 
bonds is uniform. This approximation is in fact exact 
in the Fully Segregated Model (FSM) discussed in the 
following subsection. Let us denote stretches of a = 1 
bonds by X and stretches of a = r < 1 bonds by Y. 
The two parabolas in Fig. [ll] are the steady-state J vs p 
curves for the two pure reference systems all X and all Y. 
In the disordered system, since the steady-state current 
has to be spatially constant, the possible densities are 
given by the four intersections of the line J = Jq with 
the two parabolas. If the mean density is in the range 
p < 1/2 — A (or p > 1/2 + A), then the allowed densities 
for the X and Y stretches are p\,pi (or Pi, P3) respec- 
tively. The current is in fact determined by the density 
constraint (1 — f)p\^{Jo) + f P2,?,{ Jo) = P- The variation 
of density between p\ and P2 (or p^ and p^) between X 
and Y stretches corresponds to the 'sub-bands' seen in 
Fig. (]To|a). On a macroscopic scale, however, the system 
has uniform density. Now consider what happens when 
the mean density is brought closer to 1/2. From Fig. |TJ, 
it is evident that the current would tend to increase, and 
would eventually reach the maximum allowed value J^ax 
(which equals 1/4 in the thermodynamic limit, as argued 







below). As the density is increased further, the current 
remains constant at J^axi m accordance with the max- 
imum current principle, and the excess density is taken 
care of by converting some of the X stretches of density 
pi into ones with p^ (or vice versa if p > 1/2). This 
conversion takes place adjacent to the largest stretch of 
Y bonds, leading to two macroscopic regions of differ- 
ent mean densities - one with densities (p\,p2) for the 
X and Y stretches, and the other with (p4,p 3 ). The po- 
sition of the principal shock separating these regions is 
at the location of the largest Y stretch. In the DASEP, 
the assumption of uniform density in each stretch is not 
really true, on account of the finite length of most of 
the stretches. However, the above argument provides a 
qualitative picture towards understanding the reason for 
phase separation in the DASEP. 




density (p) 



FIG. 11. Origin of phase separation in DASEP. The two 
parabolas J = rp(l — p); r = 1, 1/2, represent the J — p plots 
for the two reference non- disordered systems. 



In a certain respect, the reason for the phase separa- 
tion in the DASEP is similar to that in the single defect 
bond model studied in |ll[] — both have local segments 
in the system where the maximum allowed current is less 
than that allowed everywhere else. In the single defect 
bond model, phase separation takes place when the cur- 
rent carried by the rest of the system, with presumed 
uniform density, is larger than the maximum current al- 
lowed through the weak bond. In the DASEP, with an 
extensive number of weak bonds, the largest stretch of 
weak bonds acts as the current-limiting segment. The 
length of this stretch increases as InL with system size, 
and in the thermodynamic limit the maximum allowed 
current in this stretch is \r - equal to the maximal cur- 
rent in a pure system with only weak bonds. 

The essential point leading to phase separation is thus 
the maximum current principle, coupled with localized 
current-limiting regions in the system. In the DASEP, 
this limit is determined by long stretches of weak bonds, 
and similar considerations should apply in related mod- 
els. Consequently, we would expect a density-segregated 
phase in disordered versions of models which, in the ab- 
sence of disorder, display a maximum in the steady-state 



current as a function of density. 

E. The Fully Segregated Model 

It is useful to define a model for which many of the 
approximations made in the previous subsection are ac- 
tually exact. To this end, we study a Fully Segregated 
Model (FSM) , which is obtained from the binary random 
model above by arranging all like bonds together. Thus, 
in this model, one has two large stretches of X and Y, of 
lengths (1 — f)L and fL respectively. For the FSM in the 
thermodynamic limit, the assumption of uniform density 
within each stretch is justified, as correlations due to the 
junctions decay with increasing separation, and may be 
neglected in the bulk pi[ . 

Steady-state MC density profiles for the FSM at three 
different fillings p < 1/2 are shown in Fig. [l^ — symme- 
try under particle-hole exchange implies analogous be- 
haviour for p > 1/2 (Appendix |a|). For low densities 
(p < Pc), the two stretches have uniform bulk densities 
p x and p y related to each other by the requirement of 
equality of the two bulk currents, 

Ar(l - px) = rp y (l - p y ) = J , (32) 

and the density constraint 

(1 - f)Px + fPy = P- (33) 

These three equations determine p x , p v and Jo uniquely 
for a given p. For / = 1/2 we obtain 

__ Ap - 1 - r ± sJ{Ap - 1 - rf + 8(1 - r)p(l - 2p) 
Py ~ 2(1=7) 

p x = 2p-p y , Jo = p x (l-p x ) (34) 

This is analogous to the macroscopically homogeneous 
state of the fully random system. 
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FIG. 12. Representative density profiles of the FSM with 
r = f = 1/2 at three different fillings: (a) p = 0.24 (dotted) 
(b) p = p~ = 0.324 (solid) and (c) p = 0.4 (dashed). The 
Y stretch (r = 1/2) is in the range i £ [2000, 6000]. The in- 
set shows the variation of the bulk densities in the X and Y 
stretches as a function of the filling p. 
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As the mean density is increased, the density of each 
stretch increases, until, at a critical density p = p~ = 
\ — \-\f\. — r (the corresponding critical density on the 
higher side is p+ = 1 — p~), p y equals 1/2 (Fig. 12). 



At this density, the current equals the maximum pos- 
sible current in Y, namely, Jo = J max = r/4. As p 
is increased further, p y and Jo remain constant at 1/2 
and r/4 respectively. The density change is adjusted by 
creating a density inhomogeneity in stretch X. The two 
densities pf gh and p l ° w are related by pf gh + p l ° w = 1 
so that the currents in the two phases are equal, i.e. 
pf 9h {\ - pf gh ) = p l ° w (l - p l ° w ) = r/4. This implies 
p high,iow = (i±^/i±J)/2 (see inset in Fig. [j^). The frac- 
tion of these phases can be determined from a lever rule 
and are given by f^ah,low _ \p- p*\/(p/£9 h - fjg" 3 ). This 
locking of the density in the Y stretch at p y — 1/2 is a 
direct consequence of the maximum current principle in- 
troduced above: any change of p y from 1/2 would reduce 
the current in Y, and hence in the full system. All the 
arguments above can be applied for p > 1/2 because of 
particle-hole symmetry. Thus for p~ < p < p+ the state 
of the segregated model is analogous to the phase sepa- 
rated regime B of the random model. The size of regime 
B in the FSM is given by 2A = p+ - p~ = % /T _r 7/2. 
It closely approximates the size of the B regime in the 
DASEP (Fig. |). 



The difference between the phase boundaries of the 
DASEP and FSM Fig. [l3| comes from the fact that the 
interspersed weak-bond stretches in the DASEP have fi- 
nite lengths, and the density in these small stretches is 
not quite equal to 1/2. Close to the phase separation, we 
anticipate the mean density in a Y-stretch of length I in 
the DASEP to be of the form p Y (l) = 1/2 ± A(r)/P(r), 
with a(r) > 0. Using the distribution of the Y stretches, 
namely Py(l) — 2~ ( , we obtain the r dependence of the 
critical density 



Px 



±A{r)Y,Pr{l)PY{l) 



(36) 



where px{r « 1) ~ r is the density in the X stretches 
near the phase transition. Comparing with the phase di- 
agram for the FSM in Fig. [l3[ the correction term A(r) 
seems to be positive for all r. Further, let us suppose that 
the current in the FSM is a lower bound to that in the 
DASEP, as suggested by Fig. |l^. Then the coexistence 
curve for the random system must be quadratic near the 
critical point (J° = 1/4, p° = 1/2) - being bounded by 
two quadratics, namely, the J — p curve for the pure sys- 
tem J = p(l — p), and the coexistence curve of the FSM 
J=l/4-(l-2p) 2 . 



F. Phase-coexistence curve 

For the FSM, as r is varied we obtain different Jo vs p 
curves. The phase-coexistence curve in the current- 
density plane in the parametric form 

Jc = \, p c = \±\^/Y~r (35) 

which is the parabola J c = 1/4 — (1 — 2p c ) 2 in Fig. |l3|. 



0.25 



J 




P 

FIG. 13. The phase coexistence curve for the FSM and the 
DASEP for / = 1/2. The solid parabola is the coexistence 
curve of the FSM. The circles and the triangles are respec- 
tively the MC and mean-held phase-coexistence curves for the 
DASEP. The dashed parabola J = p(l — p) demarcates the 
allowed region for the DASEP. 



G. Correlations in the steady state 

Figure |lj shows the Monte Carlo results for the site- 
averaged density-density and height-height correlation 
functions G(Ar) and T 2 (r) in both the homogeneous and 
density-segregated regimes of the DASEP. G(Ar) is seen 
to decay rapidly over a few lattice spacings, accounting 
for the success of the mean held approximation. It is 
found that T 2 (r) grows as r — implying a roughness ex- 
ponent a = 1/2. 
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FIG. 14. Height- height correlation function F 2 (r) for 
DASEP: (a) in Homogeneous (circles), and (b) Segre- 
gated-Density (triangles) regimes. Inset shows the site av- 
eraged density-density correlation function G(Ar) defined in 
(151). The small negative values at large |Ar| arise due to the 
finite size of the system. 
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IV. DASEP WITH BACKBENDS 



As discussed at the beginning of Section III , the intro- 
duction of randomness in the easy direction of individual 
bonds alters the properties of one-dimensional disordered 
exclusion process in a crucial way. We study this in this 
section. 

The model is defined as follows: Assign quenched ar- 
rows (pointing either right or left) independently to each 
bond of a periodic chain, with probability f < \ for left 
arrows, and 1 — / for right arrows. An arrow defines 
the easy direction of hopping on each bond: a particle- 
hole exchange across a bond occurs with rate w(l + g) 
if the particle moves along the direction of the arrow, 
and w(l — g) if it moves opposite to the arrow. Since 
/ < h, there is an overall tendency for particles to cir- 
culate rightward, and the question is whether there is a 
nonzero current even in the thermodynamic limit. 

The model represents a system of hard-core particles 
in a random potential with a downward tilt. A conglom- 
eration of left pointing arrows constitutes a backbond, 
where the potential climbs up before going down again. 
Within mean-field theory it is possible to obtain an upper 
bound Ji on the current that can be carried by mutually 
excluding particles through a backbend of length I H . To 
this end, consider biased diffusion of hard core particles 
in the segment [0, 1] of a 1-d lattice, with the 'optimal' 
boundary conditions p(0) = 1 and p(l) = 0; these bound- 
ary conditions force the largest possible current through 
the segment, opposite to the bias. The master equation 
that describes transport is invariant under interchange of 
particles and holes and simultaneous relabelling of sites 
in reverse order, i.e. rij — > ni-j The boundary conditions 
respect this symmetry, implying that the steady-state 
density p(j) at site j satisfies p(j) = 1 — p(l — j). Thus in 
the steady state the number of particles in the backbend 
is 1/2 irrespective of the strength of the bias g. The prin- 
cipal effect of increasing g is to sharpen the region which 
marks the transition from the particle-rich half of the 
backbend to the hole-rich half. The steady-state profile 
approaches a step function centered at j = 1/2 as g — > 1. 

The current in the steady state is the number of parti- 
cles crossing site I in unit time. Results of a Monte Carlo 
study H are consistent with the large-^ asymptotic be- 
haviour 



./, - cxp ( -i l/A(g) 



(37) 



where A(g) is a bias-induced length given by A~ x (g) = 
ln{(l + g)/(l — g)}. This can be seen by writing the cur- 
rent within a mean field approximation as J = W(l + 

9)Pj{ 1 - Pj+i) ~ W(l - g)Pj+i( l ~ Pj), and finding the 
value of J for which the boundary conditions po = 1, 
pi = hold. For I » A(g) » l,this leads to 
J 2g e~ la ||, m agreement with ( |37| ) when g is small. 

The origin of the factor \ in the exponent in (|37]) has 
been discussed in ||, and we recount the argument in 



brief. The transport of a single particle through the back- 
bend involves two (approximately) causally independent 
steps which occur in parallel: (i) the topmost particle 
(located at site k w 1/2 in large fields) has to be acti- 
vated a distance 1/2, which requires an activation time 
r i/2 ~ exp (^ 1/ A(g)j and, (ii) the consequent hole that 
remains in the steady-state distribution moves to the bot- 
tom and is filled up, by moving each of 1/2 particles up 
through a lattice spacing. The time required is T\/2 again. 
The current J is thus proportional to t^, 2 , and conse- 
quently follows (p^). 

Since, for fixed g, the largest current that can flow 
through a long backbend (I >> A(g)) is exponentially 
small in its length, the largest current through the 1-d 
lattice of length L is determined by the length l*(L) of 
the largest backbend encountered. Since the probabil- 
ity of occurrences of I consecutive left-pointing arrows 
on bonds is proportional to / , we may estimate I* from 
Lf l = C where C is a constant of order unity. Substitut- 
ing in ( |37|) , we find that the current falls with increasing 
lattice size as 



J(L) 



-he 



(38) 



with O^ 1 = A(g)ln/. Thus the current is expected to 
decay as a power law in L, with a bias-dependent power, 
and to vanish in the thermodynamic limit. Figure |l5| 
shows the result of Monte Carlo simulations. 
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FIG. 15. Size dependence of the steady-state current in the 
backbend model (DASEP with some reversed bonds) for two 
sets of parameters: (a) g = 0.33, / = 0.25 (triangles) and (b) 
r = 0.54, / = 0.3 (circles). Each point represents an average 
over 40 realizations of disorder in (a) and 100 realizations in 
(b). The straight lines have slopes of —9/2 with 9 — 0.5 and 
1.0 respectively as predicted by (|3^). 

As with the milder sort of disorder discussed in Sec- 
tion 111, the state is strongly inhomogeneous and shows 
macroscopic regions of high and low density. Figure |l]c 
shows the time-averaged density profile for a typical con- 
figuration of bonds. There is a large shock around the 
rate-limiting backbend, which separates the two regions. 
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For fixed lattice size L with an associated longest back- 
bend l*(L), the current is a nonmonotonic function of g. 
This can be seen as follows. If g is low enough that 
A(<?) >> l*(L), linear response theory would imply that 
current J grows linearly with g. On the other hand, if g 
is large enough that A(g) << l*(L), the current falls with 
increasing g according to (|37]). In between, J achieves a 
local maximum when A(g = g ma x) — which im- 

plies g max ~ 1/ln-L. 

The argument given above implies that the current 
carried by a system of hard-core particles through the 
randomly backbending lattice vanishes in the thermody- 
namic limit, no matter how small the bias g. This is in 
contrast to the behaviour of noninteracting particles in 
the same random environment, where the drift velocity 
vanishes only if the bias is strong enough J2^,||. The 
difference can be traced to the possibility, in the non- 
interacting case, of compensation by a large build-up of 
density at the bottom of a backbend, which then succeeds 
in driving a finite current over the backbend. This option 
does not exist once repulsive hard-core interactions come 
into play, and the current vanishes in the thermodynamic 
limit. 



V. CONTINUUM DESCRIPTION 

It is interesting to ask whether the behaviours found 
above in the disordered lattice gases can be reproduced 
using a continuum description of the problem. Though 
we have not pursued this question to its logical end, we 
discuss in this section some general constraints that a 
continuum description should satisfy. 

The steady-state of the disordered system of interact- 
ing particles is described by a spatially varying time- 
averaged density profile po(r) = (n(r)). The time evolu- 
tions of fluctuations around the mean density profile are 
governed by the continuity equation 

J^p(r,t)+V- J(p,r,t) = 0. (39) 

The coarse-grained current J may be phcnomcnologically 
divided into three parts 

J{r, t) = J sys {p{r, t),r) + J diff + e{r, t) . (40) 

A local hydrodynamic assumption has been made in writ- 
ing the systematic part of the current, J sys as a function 
of the local density, p(r,t). The explicit r dependence of 
J sys comes from the breaking of translational invariance 
by quenched disorder, and its exact form of J sys depends 
on the microscopic dynamics of the model. The diffu- 
sive current, Jdiff = D(r) ■ V[p(r,i) — Po( r )L involves 
a quenched diffusion tensor, D(r). The noise, e(r,t), is 
to mimic, on a mesoscopic scale, the stochastic nature 
of the evolution. It is usually assumed to be Gaussian 



distributed and 5-correlated in space and time with van- 
ishing spatial and temporal averages. 

To obtain the time evolution of the density fluctua- 
tions p = p(r,t) — po(r), we expand J sys in powers of p 
as J sy s{po(r),r) + c{r)p + X(r)p 2 ■ ■ ■ and put it in (|3S|). 
This results in 

d t p = V • [D(r) ■ Vp - c(r)p - \{r)p 2 e(r, t)] (41) 

In one spatial dimension, the above reduces to (with r 
replaced with x) the form 

dtp = d x [D(x)d x p - c(x)p - X(x)p 2 ■ ■ ■ - n{x, t)} (42) 

which was considered in In steady state, the time- 
averaged current must be independent of a;. As both 
Jdiff and e(x, t) vanish under time averaging, the con- 
straint to be satisfied is 

d x (J sy s(x)) = (43) 

which is important to account for, as the coefficients in 
J S ys are explicitly space-dependent, i.e., 9 x [A(a;)(/5 2 ) + ...] 
must vanish. 

In their attempt to study a continuum model which 
describes the DASEP, Becker and Janssen (BJ) || write 
the current J' in powers of 4>(x, t) = p(x, t) — p, the den- 
sity fluctuation away from the overall particle density p 
in the system. In one dimension, the form quoted || for 
the DASEP is 

J' = (1 - 2p)4>(x, t) - 4> 2 {x, t) + n(x) (44) 

where rj(x) is an additive quenched noise term. Since 
the time-averaged density po(x) varies in space, it is ev- 
ident that 4> has a nonzero expectation value (<j)(x, t)) = 
Po(x) — p. As we have seen in Section III, po(x), and 
thus ((j>(x,t)}, can show strong variations, especially in 
the density-segregated regime. Spatial constancy of the 
average current demands that the time average of the 
right hand side of (Q) must satisfy 

cU(l - 2p){4>{x, t)) - (c)> 2 (x, t)) + r)(x)] = 0. (45) 

Even when this condition is satisfied, it is not entirely 
clear that this continuum theory actually represents the 
DASEP in one dimension. BJ argue that if dJ/dp ^ 0, 
then disorder is irrelevant on the large scale, a conclu- 
sion that is supported by [p|j26|. However, their conclu- 
sion, that the condition dJ/dp = holds only if p = 4, 
does not seem to be correct, at least for the DASEP in 
one dimension; we have seen in Section IV, there is an 
entire range of densities (5 — A < p < | + A) where 
dJ/dp vanishes, associated with segregation of density 
on a macroscopic scale. 
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VI. EQUIVALENT INTERFACE MODELS IN 
ONE DIMENSION 



In 1-d, both the DASEP (with or without backbends) 
and the GDDP - for which the maximum occupancy per 
site is 1 - are equivalent to stochastic growth models of 
a 1-d interface moving in a 1-d medium. Corresponding 
to each particle- hole configuration {rii} is assigned an in- 
terface profile {Hi} through H t = £^(1 ~ & 
Pictorially this means that each particle corresponds to 
a —45° downward line segment, while a hole corresponds 
to an upward one (Fig. [l6]). Thus, clusters of particles 
and holes translate to ±45° slope segments, and the inter- 
face has a mean slope which vanishes when the particle 
density is 1/2. Away from half-filling, periodic bound- 
ary condition for the lattice becomes a helical boundary 
condition for the interface. Junctions between adjacent 
particle and hole clusters correspond to corners in the 
interface profile. 



along the interface - being proportional to the spatially 
constant steady-state current of the corresponding parti- 
cle model. 

We now discuss the various qualitatively different 
regimes which arise in the interface growth models. Fig- 
ure 17 shows time averaged steady-state interface profiles 
(Hi) corresponding to the three regimes of driven parti- 
cle systems illustrated in Fig. [I] in Section ffl. In all the 
three cases we started from an initially uniform profile 
- corresponding to a random distribution of particles on 
the lattice. 





V 


a 









b 





L 



H. 




(a) Toom interface Model 




(b) Single-step Model 



FIG. 16. Mapping between driven particle systems in 
d — 1 and growing interfaces, (a) Toom Interface dynam- 
ics corresponding to the GDPP, (b) Single-step model corre- 
sponding to the DASEP. 

Evolution of the interface is dictated by the dynamics 
of the corresponding particle system. The GDDP corre- 
sponds to the slice-wise motion of segments of a Toom 
interface in the low noise limit p0| , while the DASEP 
corresponds to the corner-flip 'single step' growth model 
|Q (Fig. |l6|). In both cases particle movement to the 
right (or hole move to left in GDDP) corresponds to lo- 
cal forward growth (deposition) of the interface while 
a leftward particle move (e.g. in DASEP with back- 
bends) corresponds to local backward growth (evapora- 
tion). The quenched jump rates for the particle moves 
implies columnar disordered growth-rates for the inter- 
face: the normal growth rate of the interface at a fixed i 
is the same irrespective of the location of the interface at 
successive times. In the long time limit, the mean local 
forward speed of the interface is the same at all points 



A 




FIG. 17. Interface morphology in the three regimes: (a) 
Uniformly moving interface with a uniform slope, (b) inter- 
face with large segments of different slopes moving with a 
non-vanishing speed, and (c) interface with large segments of 
different slopes moving with a speed which vanishes in the 
thermodynamic limit. The insets show blow-ups of regions 
enclosed in the dashed boxes. 



Figure |l7| a is an interface with a non-zero mean tilt, 
which has net forward growth rate at all points. The in- 
terface has a uniform slope on a macroscopic scale and 
moves with a finite non-zero speed preserving its mean 
tilt. This case corresponds to the Homogeneous regime 
depicted in Fig. ffla. On a microscopic scale the interface 
has frozen-in roughness (Fig. |l7|a, inset) corresponding to 
the microscopic shocks in the steady-state density profile 
of Fig. |a. 

If the mean tilt vanishes, but the interface still has 
a net forward growth rate at all points, then the initial 
uniform profile at t — coarsens into large segments of 
different mean slopes at long times (Fig. |l7jb). These 
segments have frozen roughness on microscopic scales, 
similar to the non-zero tilt case (Fig. 0b, inset). The 
interface moves with a finite speed preserving its mean 
shape and mean vanishing tilt. This corresponds to the 
Segregated-Density regime of Fig. [ljb. 

In addition to the frozen roughness on the microscopic 
scales, we can define the kinetic roughness as the equal- 
time mean-square height fluctuations around the steady- 
state profile. We consider the zero-mean height variables 
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hi(t) = Hi(t)—(Hi) defined in (|l7|), and define the rough- 
ness exponent a through ((hi +r — hi) 2 ) ~ r 2a . As dis- 
cussed in Subsections ||D,E and III F, the roughness ex- 
ponent a = 1/2, in both the cases above. 

In Fig. |l7|c, which corresponds to the Vanishing- 
Current regime of Fig. [l]c, the profile resembles that cor- 
responding to the Segregated-Density regime in that it 
has large segments of different mean slopes. However, 
in this case the interface is stationary in the thermody- 
namic limit - reflecting the vanishing of the steady-state 
current in the particle system. In interface language, this 
situation can be visualized as a case of interface growth 
where there are local stretches of the interface having 
a net backward growth (evaporation) rate. Though, on 
the average, there are more forward-growth regions than 
backward-growth ones, in the limit of large system size, 
arbitrarily long evaporation stretches effectively pin the 
interface. 

Turning now to continuum description of the dynam- 
ics of these interfaces, at least in the cases where the 
mean speed of growth is non-zero, the sum in the def- 
inition of the height variable is replaced by an integral 
of the coarse-grained particle density, H(x,t) = J [1 — 
2p(x', t)]dx' . The growth equation for H (x, t) is governed 
by (spatially) integrated one dimensional version of the 
continuity equation (|3|): d t H(x,t) = J(d x H,x,t). The 
fluctuations h{x,t) = H{x,t) — Hq(x)\h H around the 
steady-state profile H (x) = (H(x, t)), are governed by 

d t h = D(x)d xx h - c{x)d x h + ^\(x)d x h 2 h 2n(x, t) 

(46) 

obtained by integrating (^) . The absence of any additive 
quenched spatial noise term in ( fl6|) is due to the spatial 
constancy of the growth speed of the interface dictated 
by the same constraint on the steady-state current. In 
this respect fl46| ) differs from the model discussed in |27]] 
where such a term arises naturally due to absence of any 
such constraints. As can be readily verified by power 
counting, an additive quenched columnar term is highly 
relevant in the RG sense and leads to much rougher in- 
terfaces than the a = 1/2 interfaces described by (fig). 



VII. CONCLUSION 

In this paper, we have studied the stationary current- 
carrying states of driven lattice gas models with quenched 
disordered hopping rates. The principal results are of two 
types — first the exact determination of the steady states 
for a class of disordered models, and second the demar- 
cation of distinct regimes of behaviour on macroscopic 
length scales, as a result of disorder. In this section, 
we briefly review these results, and discuss some related 
open problems. 

The steady states of a family of disordered models - 
the disordered drop-push process (DDPP), and related 



models — have been found in all dimensions by an ap- 
plication of the condition of pairwise balance. The result 
is a product measure state, with site-dependent weights, 
reflecting the microscopic disorder in the model. The cur- 
rent has been computed as well. The system is character- 
ized by a strictly uniform current density, and a coarse- 
gained particle density that is approximately uniform. 
On a macroscopic scale, the state is homogeneous. 

Disorder can lead to macroscopically non-homogeneous 
states, as in the 1-d disordered asymmetric simple ex- 
clusion process (DASEP). Our numerical and mean-field 
results show that a macroscopically density-segregated 
state occurs in the DASEP model with no backbends, 
for densities in a finite range around half-filling. The 
origin of density separation is traceable to the existence 
of a largest current that can be carried by a stretch of 
weak bonds. This low current can be sustained in the 
rest of the lattice only by separating the density into dis- 
tinct large and small values in macroscopic regions of the 
lattice. 

Backbends introduce a third type of possible behaviour 
in one dimension. Like the stretch of weak forward bonds 
in the density-segregated regime, a backbend rate-limits 
the current, leading to density segregation. However, 
there is an important difference: the longer the stretch of 
weak forward bonds, the closer is the current to a finite 
asymptotic value; by contrast, the longer the stretch of 
reverse-biased bonds, the smaller the current — it de- 
creases exponentially fast with backbend length. Since 
the probability of occurrence of a backbend decreases ex- 
ponentially with its length, the result is a current that 
decreases as a bias- dependent power of the overall size of 
the lattice. 

The crucial physical point which underlies the be- 
haviour in each of these regimes is the requirement that 
the steady-state current be constant at all points in the 
system. Continuum field-theoretical approaches must en- 
sure that this local constraint is respected; while this is 
automatically assured for translationally invariant sys- 
tems, it may need special care to guarantee in disordered 
systems. 

It would be of interest to generalize our results to 
higher dimensions and also to include interactions be- 
tween particles at different sites. A few scattered results 
along these lines are available. 

(i) For the drop-push class of problems, we have seen 
in Section || that the exact steady state even in higher 
dimensions is characterized by inhomogeneous product 
measure. On large scales, this leads to a homogeneous 
state. 

(ii) The transport of particles with hard-core interactions 
through the infinite cluster of a randomly diluted lat- 
tice above the percolation concentration has been studied 
HI). In a certain regime of dilution, backbends act as lo- 
cal traps, but unlike the one-dimensional case considered 
in Section IV, there exist infinitely long paths on which 



the length of every backbend is less than a fixed value 
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The sub-network of such paths is expected to 
carry a current which then remains finite in the thermo- 
dynamic limit. There is thus no vanishing-current regime 
in this system. 

(iii) With attractive interactions between particles, the 
driven lattice gas system with nearest-neighbour hopping 
is known to undergo phase separation below a certain 
temperature [p9pf| . A numerical simulation showed that 
the addition of a low concentration of blocked sites did 
not alter the critical behaviour of this system JtJ . 

More systematic studies of higher-dimensional systems 
are called for. In particular, it would be interesting 
to know whether disorder-induced large-scale inhomo- 
geneities, akin to the phase separation found in one di- 
mension, persist in higher dimensions as well. 
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APPENDIX A: SYMMETRIES OF THE 
CURRENT IN THE DASEP 

In this appendix we discuss two types of symmetry 
transformations in the DASEP which leave the steady- 
state current invariant, (i) The first, which involves flip- 
ping the directions of all jumps and interchanging all par- 
ticles and holes, can be proved exactly. Under these oper- 
ations, the steady-state weights of particle configurations 
and the magnitude of the current are shown to remain 
invariant, (ii) The second type of symmetry is restricted 
to l-d systems in which only one way jumps are allowed 
on each bond. The symmetry transformation consists of 
flipping the direction of jump on each bond. In this case, 
we have convincing numerical evidence that the current 
is invariant, though there appears to be no simple rela- 
tionship between steady-state weights of various particle 
configurations. We have no general proof of this result. 

(i) Easy direction flip and particle-hole interchange: 
Consider the DASEP, with quenched disordered, unequal 
forward and backward hopping rates on each bond. Be- 
low we explicitly deal with the one-dimensional case, but 
the results are easily generalizable to higher dimensions. 
Let R denote a particular realization of disorder and R 
the realization obtained from R by flipping the easy di- 
rection of jumps on all the bonds: If R denotes the set 
of transition rates {V^ 1 -' }, then R corresponds to the set 

{W* 3 = W° 1 }. Also, let us denote by C and C two par- 
ticle configurations related to each other by particle-hole 
interchange. Clearly, if C is an allowed configuration of 
the system at filling p, then C corresponds to a filling 
1 — p, and thus there is a one-to-one correspondence be- 
tween the configurations at the two fillings. Now, let 
be the configuration obtained from C by exchanging the 
occupation numbers at the sites i and j. It is easy to 
see that the two transition rates, W P (C — » C 4 - 7 ) in real- 
ization R, and W\- P {C — > C y ) in realization R are equal 



(Fig. |l8|), i.e. the transition matrices W p and Wi- P in 
the two realizations have identical entries. 
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FIG. 18. Invariance of the current under R —* R and 
p — > 1 — p for the DASEP in l-d. Jump directions on ev- 
ery bond are reversed from R to R. Configurations C,C are 
related by particle-hole interchange. So are C %] and C iJ ' . 

Identification of the two W-matrices implies that the 
invariant measures P p and P\- p satisfy 



P P (C) = Pi-p(C). 



(Al) 



Using this, together with the identity rij(C) = 1 — nj(C), 
we can relate n-point correlations at fillings p and 1 — p. 
In particular, the site densities at the two fillings are re- 
lated as 



i)p,R = 1 — 



'l-p,R- 



(A2) 



Further, the steady-state currents in R and R at the two 
fillings are equal in magnitude: 

J (p, R) = W^( ni (l - n,))p,fi - W ji ( nj (l ~ ni )) PtR 



-J (l-p, R). 



(A3) 



The negative sign in the last step reflects the fact that the 
direction of the current is opposite in the tw o systems. 

In the case of the FSM (subsection [II E ) the realiza- 
tions R and R are identical (apart from a reflection). 
Hence in this case the above arguments imply that the 
magnitudes of the currents at the two fillings p and 1 — p 
are equal. A similar result holds for the single defect 
bond case studied by Janowsky and Lebowitz 



For the DASEP, with or without backbends, (A3) has 
the corollary that the disorder averaged currents at the 
two fillings p and 1 — p are equal. 

(ii) Easy direction flip only (for l-d, forward hopping): 
Above, we found that the current is invariant when the 
filling is changed from p to 1 — p, provided that the disor- 
der realization is changed from R to R. Here we observe 
(based on numerical evidence) that in l-d, with only for- 
ward hopping, the result is also true for realization R on 
its own, i.e. 



Jo(p,R) = J (l-p,R). 
In view of (|A^) this is equivalent to 
Jo(p,R) = -Mp,R). 



(A4) 



(A5) 
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The claim is easily verified for the single particle (or 
single hole) case using an explic it fo rm for the current: 

Jo = (E^m+iT 1 m- Also ' © is true at P = V 2 
since (A4) is an identity at this filling. 



We do not have a proof for (A4) or (A5) in the general 
case, but they seem to be borne o ut n umerically. For 
instance, we studied the validity of (A5) for a system of 
size L = 6 with N = 2 particles. We took R to be the set 
= 1,...,6} = The invariance 

of the current is verified up to 1 part in 10 8 . We also 
studied the steady-state probabilities of all the 6 C2 par- 
ticle configurations for each of R and R, both by Monte 
Carlo and by Lanczos iteration of the stochastic evolution 
operator. There seems to be no straightforward corre- 
spondence between the two sets of steady-state weights. 
This suggests that there should be a proof of the invari- 
ance of the current which does not rely on identifying the 
weights of individual configurations. 
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